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Abstract 

In this work, we have done a completely microscopic calculation using a many-body variational 
method based on the cluster expansion of energy to compute the asymmetry energy of nuclear mat¬ 
ter. In our calculations, we have employed the AVi^ nuclear potential. We have also investigated 
the temperature and density dependence of asymmetry energy. Our results show that the asymme¬ 
try energy of nuclear matter depends on both density and temperature. We have also studied the 
effects of different terms in the asymmetry energy of nuclear matter. These investigations indicate 
that at different densities and temperatures, the contribution of parabolic term is very substantial 
with respect to the other terms. Therefore, we can conclude that the parabolic approximation is a 
relatively good estimation, and our calculated binding energy of asymmetric nuclear matter is in 
a relatively good agreement with that of semi-empirical mass formula. However, for the accurate 
calculations, it is better to consider the effects of other terms. 
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I. INTRODUCTION 


Nuclear matter is an extremely large hypothetical system of constrained nucleons which 
interact through the strong nuclear force, and can be considered as an ideal model for 
the nuclear matter within the heavy nuclei |^. The analysis of heavy nucleus behavior 
in high density and hnite temperature is a challenging issue in modern nuclear physics 
[^. The nuclear matter asymmetry energy is an important subject in nuclear physics and 
astrophysics. According to research conducted js-lO], the nuclear matter asymmetry energy 
plays a crucial role in determining many important nuclear properties such as the neutron 
skin of nuclear system, structure of nuclei near the drip line, and neutron star mass and 


radius. They are especial important constraints on_t 
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le nuclear matter equation of state and 


nl,Q. 


The nuclear matter equation 


the density dependence of the asymmetry energy 
of state plays a central roles in the study of high-energy in nuclear physics, especially in 
heavy ion collisions. The nuclear matter equation of state at zero temperature governs the 
structure of cold neutron stars, the equation of state for hnite temperatures is necessary for 
ies of core collapse supernovae, black hole formation and cooling of the neutron stars 


stu( 
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Due to the strong interaction between the nucleons, we must completely include the 
nucleon-nucleon potential in the nuclear matter calculations. The nuclear potential is usu¬ 
ally made out by htting the nucleon-nucleon scattering data and deuteron properties. The 
potential models which obtained only from htting the np scattering data frequently oher 
a weak description of pp and nn data [^. Recently, a nuclear potential has been intro¬ 
duced by Wiringa et al. which extracted by htting np and pp data as well as low energy nn 


scattering parameters and deuteron properties 


17l | . This potential is written in an operator 


format which depends on the value of J, L, S and T as well as of nucleon-nucleon pair. 
Therefore, in the calculations with this potential, a completely microscopic calculation in 
which the third component of isospin (T^) would be explicitly considered is required. This 
is a very important point in the calculations for asymmetric nuclear matter consisting of Z 
protons and N neutrons, where in general Z ^ N. The microscopic studies on asymmetric 
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nuclear matter are very limited, and usually the asymmetry energy of nuclear matter is 
computed employing a parabolic approximation. Therefore, it is interesting to study the 
validity of this estimation for the asymmetric nuclear matter. 

One of the powerful microscopic methods in the nuclear matter calculations is the vari¬ 
ational method based on the cluster expansion of energy. In recent years, we have done 
several works on the nucleonic matter using this method 
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25| . The purpose of present 


work is to compute the asymmetry energy of nuclear matter using AVis potential at hnite 
temperature employing the mentioned microscopic technique. We also verify the validity 
of semi-empirical mass formula at different temperatures and densities. The dependence of 
asymmetry energy of asymmetric nuclear matter on density and temperature is also inves¬ 
tigated. 


II. FORMALISM 


Here, a system with A interacting nucleons (including N neutrons and Z protons) is 
considered in order to calculate the binding energy of asymmetric nuclear matter with AVis 
potential. As it was mentioned in the preceding section, we use the variational method 
based on the cluster expansion of the energy. In this method, at hrst we choose a trial wave 
function as follows, 

ij = F<h, (1) 

where $ is the antisymmetric wave function of A non-interacting nucleons and F = 
F{1... A) is the inter-nucleon correlation function included by the interaction between the 
nucleons. We use the Jastrow ansatz for the correlation function. 


i>j 


( 2 ) 


where S is the symmetrizing operator and f{ij) is the two-body correlation function. Using 
the above trial wave function included by this correlation function, we get a cluster expansion 
for the energy which up to the two-body term is as follows. 


£■(!/]) 


A <^|^> - 


( 3 ) 


At hnite temperature, the one-body term Fi for an asymmetric nuclear matter consisting 


3 






of Z protons and N nentrons is as follows, 


El = + e[''\ 


(4) 


Labels p and n are nsed for the protons and nentrons, and is. 


2mh) 

where n^'^\k,T, is fermi-Dirac distribntion fnnction, 

1 


(5) 




( 6 ) 


2m(') 


Exp{l3[e^‘^\k) — -|- 1 

In the above eqnation, /i^*^ is the chemical potential, jg the nnmber density and = 
is the single particle energy of nncleon z, and /3 = -j^ (T is the temperatnre). The total 
density of system is, 

p = ph)+pW. (7) 

For the asymmetric nnclear matter, we can also dehne the asymmetry parameter as follows, 

N - Z p^^'> - 


^ N + Z p^^) + p(p) ■ 

This dehnition leads to the following relations for the nncleon density. 


( 8 ) 


= (1 - Op, p^"^ = (1 + Op- 


(9) 


By comparing Eqs. (jD]), (jS]) and (jl]), we can conclnde that the energy of asymmetric nnclear 
matter is different for varions asymmetry parameter. 

The two-body energy has the following form, 

^2 = ^ I ^^(12) I ij - ji >, (10) 

ij 

where 

>’(12) =-^1/(12), |V,^/(12)|] + /(12)V'(12)/(12). (11) 

B(12) is the nncleon-nncleon pair potential. The two-body correlation fnnction, /(12), is 
considered as follows, 

/(12) = ^/l'’>(ri2)0"’>(12), (12) 

p=l 
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where operators 0 ^\l2) are, 


Qip 1-3)('^2) = (i-iai.cJa), (^ + + ^^ 12 ), 

(4 + Y^ai.a2--^i2). 


(13) 


In above equation, S '12 = 3((Ti.f)(cr2.r) — ai.a2 is the tensor operator. We employ the AVis 
nuclear potential in our calculations. This potential has the following form 17l |. 
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1 /( 12 ) = 

p=i 


(14) 


The first fourteen operators are. 




= 1-14 


12 


= l,ai.(T2,ri.r2, (o-i.(T2)('ri.r2), S'12, S'i2(ri.r2), L.S, 

L.S(ti. r 2 ), L^((Ti.0-2), (ti. r 2 ), L^((Ti. 0-2) (ri.r 2 ), 
(L.S)^(L.S)2(rl.r2), 


(15) 


and four additional operators which break the charge independence of nuclear force are. 




= 15-18 


12 


— T12, Ti2(cri.cr2), 5 'i 2 Ti 2 , (Tiz + T2z)- 


(16) 


Ti 2 = 3(ri.r)(r2.r) — T1.T2 is the isotensor operator. 

Using the above formalism, we can compute the binding energy of asymmetric nuclear 
matter. The procedure of these calculations has been fully presented in our previous articles 
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III. RESULTS AND DISCUSSIONS 

We have calculated the binding energy of asymmetric nuclear matter using the method 
discussed in the previous section at various densities and temperatures. Our results for 
the binding energy versus asymmetry parameter (^ = j^~^ ) have been shown in Fig. [1] 
for different densities and temperatures. For all densities and temperatures, it is seen that 
the binding energy increases as the asymmetry parameter increases. We have found that 
for each value of temperature and asymmetry parameter, below a specific value of density, 
the binding energy decreases by increasing density, while above this density, it increases as 
density increases. This specihc density is called the saturation density which depends on 
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both temperature and asymmetry parameter. For each density, by comparison of the energy 
curves, we observe that by increasing the temperature, the slope of energy curve versus 
asymmetry parameter decreases. Fig. [T]shows that for each temperature, the increasing rate 
of binding energy versus asymmetry parameter has different values for different densities. 

For each value of density at a specihc temperature, we can write our results for the binding 
energy versus asymmetry parameter ((^) as a polynomial function in term of different powers 
of By comparing the extracted polynomial function with the semi-empirical mass formula 
used for the binding energy of asymmetric nuclear matter {E = ay + agym^'^), the requirement 
of presence of other powers of ^ can be analyzed, and it is determined whether the semi- 
empirical mass formula would be broken or not. We have found that the best function which 
is htted with our results for the binding energy of asymmetric nuclear matter is a polynomial 
up to the fourth power of 

E = Uq 0 - 2 ^^ “1“ 

The coefficients of various powers of ^ (ao, ai, 02 , 03 and 04 ) versus density at different 
temperatures have been given in Tables m and mil By investigating the values of Oq from 
these tables, we observe that for each temperature, ao decreases by increasing the density up 
to a specihc density, and then above this density, Oq increases. By comparing these tables, 
we can hnd that for each density, with an increase in temperature, an increasing in oq is seen. 
From Tables min] and mil by investigating the coefficients Oi, 03 and 04 , we can conclude 
that these coefficient are density and temperature dependent. The given values of 02 in 
these tables indicate that for each temperature, the value of 02 increases by increasing the 
density. We have found that at high densities, the increasing rate of 02 is lower than that 
at low densities. A comparison between the values of 02 in different temperatures shows 
that for all densities, the value of 02 decreases as the temperature increases. Here, we can 
conclude that the value of 02 depends on both density and temperature. 

Now, by more investigating Tables m El and IIIIl we see that for all densities and tem¬ 
peratures, with respect to the other coefficients, the coefficients Oq and 02 have the main 
contributions in the binding energy. However, we test the the effects of ai, and 04 relative 
to a 2 as follows. By obtaining the magnitude of ratio it is clear that at low densities, 
this ratio is of order 10 “^, therefore the effect of Oi in the binding energy of asymmetric 
nuclear matter is negligible. We have found that at high densities, this ratio is of order 10”^, 
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TABLE I: ao, ai, 02, 03 and 04 (MeV) versus density (p) for temperature T = 0 MeV. 


p 

ao 

ai 

02 

as 

04 

0.05 

4.921 

-0.05167 

12.01 

-0.117 

2.574 

0.1 

-10.08 

-0.01954 

17.64 

1.548 

0.4468 

0.15 

-13.76 

-0.1445 

23.76 

0.1695 

1.371 

0.2 

-16.3 

-0.2257 

28.82 

-0.6103 

2.067 

0.25 

-17.85 

-0.3245 

33.57 

-1.974 

3.028 

0.3 

-18.45 

-0.4127 

37.37 

-1.949 

3.02 

0.35 

-18.19 

-0.4084 

40.68 

-1.887 

3.007 

0.4 

-17.1 

-0.5143 

44.33 

-2.762 

3.524 

0.45 

-15.23 

-0.6283 

47.94 

-3.658 

3.932 

0.5 

-12.61 

-0.6323 

50.63 

-2.82 

3.355 


TABLE 11 : oq, ai, 02, 03 and 04 (MeV) versus density (p) for temperature T = 10 MeV. 


p 

ao 

ai 

02 

as 

04 

0.05 

6.73012 

-0.04794 

8.94846 

0.48205 

0.75552 

0.1 

-0.02839 

-0.17483 

15.43355 

-0.33542 

0.77439 

0.15 

-4.86116 

-0.13702 

20.64022 

0.60392 

0.2858 

0.2 

-8.2671 

-0.26539 

26.04789 

-0.0697 

0.79922 

0.25 

-10.48181 

-0.24034 

30.19788 

0.29485 

0.71646 

0.3 

-11.63058 

-0.40422 

34.81898 

-0.85014 

1.4016 

0.35 

-11.79859 

-0.62335 

39.13478 

-1.82609 

1.88743 

0.4 

-11.07976 

-0.47165 

41.84846 

-0.96744 

1.5085 

0.45 

-9.50737 

-0.5509 

45.37664 

-1.40054 

1.66413 

0.5 

-7.1490 

-0.52195 

47.92776 

0.1053 

0.63481 


therefore it is better to consider the effect of ai for the accurate calculations. By comparing 
the magnitude of 03 and a 4 , it is seen that for all temperatures and densities, the magnitude 
of 04 is always greater than that of 03 . Our results show that the highest order of magnitude 
of ratio — reaches about 10“^. For the magnitude of —, that is of order 10“^. The above 
results indicate that the dominant term in the asymmetry energy of nuclear matter is related 
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TABLE III: oq, oi, 02 , 03 and 04 (MeV) versus density {p) for temperature T = 20 MeV. 


P 

ao 

ai 

a2 

as 

04 

0.05 

22.50472 

-0.21782 

8.02506 

-2.70781 

2.36938 

0.1 

16.14381 

-0.12781 

11.66631 

0.04014 

0.33926 

0.15 

11.29118 

-0.09781 

15.88062 

1.21584 

-0.27871 

0.2 

7.65327 

-0.18837 

20.78746 

0.18392 

0.49755 

0.25 

5.12343 

-0.24998 

24.86982 

0.36818 

0.36781 

0.3 

3.61406 

-0.39247 

29.43257 

-1.14565 

1.30578 

0.35 

3.07363 

-0.34095 

32.36919 

0.23839 

0.38758 

0.4 

3.45085 

-0.56422 

36.75817 

-1.36493 

1.19642 

0.45 

4.68426 

-0.55863 

39.84991 

-0.73474 

0.63449 

0.5 

6.71498 

-0.48627 

42.39283 

1.09249 

-0.69099 


to 02 , and the other terms which are related to oi, 03 and 04 are small with respect to that 
of 02 - However, in order to do more accnrate calcnlation, these terms could be considered. 

Here, for the binding energy of asymmetric nuclear matter, we rewrite Eq. mill as the 
following relation 

E Oj^ T Engyirji, (18) 

where o^ = oq is the volume energy and Easym is the asymmetry energy (for the symmetric 
nuclear matter in which N = Z, Easym is zero). 


■‘-^asym ■‘~^asym ' '‘~^asym ' ^asym ' '‘~^asym' 


(19) 


In fact the contribution of asymmetry energy is considered as four terms. Now, at a specific 
density and temperature, we compare the effects of different terms of Easym with respect to 
E^sl/m = 0-1) high = 0.9) values of asymmetry parameter. Our results 


indicate that for the low values of asymmetry parameter, the magnitude of ratio — 

^asym 

is signihcant, especially at higher densities, therefore, for this case of asymmetric nuclear 
matter, it is better to consider the effects of E^^y^. However, by increasing the asymmetry 

^( 3 ) 

parameter, the effect of this term decreases. The magnitude of ratio —and at low 

^asym ^asym 

value of asymmetry parameter, is very small. The magnitude of these ratios at high values 
of asymmetry parameter becomes more important, but they are yet small. 
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IV. SUMMARY AND CONCLUSIONS 


We have considered the asymmetric nuclear matter including N neutrons and Z protons 
which interact through the nuclear force. For this system, the binding energy has been 
computed by a fully microscopic variational many-body technique based on the cluster ex¬ 
pansion of the energy. In our formalism, we have presented the AVis two-nucleon potential. 
We have also calculated the asymmetry energy of nuclear matter. It is found that the 
asymmetry energy is a function of density and temperature. For different densities and tem¬ 
peratures, our results for the binding energy versus asymmetry parameter have been htted 
as a polynomial function. It was a polynomial of forth power of asymmetry parameter. We 
have seen that all the coefficients of this polynomial are density and temperature depen¬ 
dent. By comparing our polynomial function with the semi-empirical mass formula used 
for the binding energy of asymmetric nuclear matter, the requirement of presence of other 
powers of asymmetry parameter has been analyzed, and the validity of semi-empirical mass 
formula has been investigated. According to our results, the square term in the asymmetry 
energy of nuclear matter is the dominant term, and the other terms compared to this term 
are small. Therefore, we can ignore the effects of non-square terms. Our study shows that 
the asymmetry energy of asymmetric nuclear matter has a relatively good agreement with 
the semi-empirical mass formula. However, according to our calculations, to increase the 
accuracy of calculations, it is better to include the other terms in the asymmetry energy. 
For example, the hrst order term at low values of asymmetry parameter is considerable, 
especially at low densities. 
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FIG. 1: The binding energy per nucleon of asymmetric nuclear matter versus asymmetry parameter 
(^) for densities p = 0.05 (full curve), 0.25 (dashed curve) and 0.5 fm~^ (dotted curve) at different 
temperatures (T). 
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